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We present results from our analysis of Chandra X-ray Observatory, W. 
M. Keck Observatory, and Karl G. Jansky Very Large Array (VLA) images of 
the Crab Nebula that were contemporaneous with the 7-ray flare of 2011 April. 
Despite hints in the X-ray data, we find no evidence for statistically significant 
variations that pinpoint the specific location of the flares within the Nebula. The 
Keck observations extend this conclusion to the "inner knot", i.e., the feature 
within an arcsecond of the pulsar. The VLA observations support this conclu- 
sion. We also discuss theoretical implications of the 7-ray flares and suggest that 
the most dramatic 7-ray flares are due to radiation-reaction-limited synchrotron 
emission associated with sudden, dissipative changes in the current system sus- 
tained by the central pulsar. 



1. Introduction 

The Crab Nebula, the relic of a stellar explosion recorded by Chinese astronomers in 
1054, has a special place in the history of astronomy. It is our most frequently observed 
laboratory for high-energy astrophysics. Located at a distance of ~ 2 kpc, the system is 
energized by a pulsar of spindown luminosity L v \ sr « 5 x 10 38 erg/ s and current spin period 
P « 34 ms. The history and general properties of the system are nicely summarized in the 
review by Hester (2008). Optical and X-ray images (Hester et al. 1995; Weisskopf et al. 2000; 
Hester et al. 2002) of the inner nebula show features such as an inner ring, toroidal structure, 
knots, and two opposing jets originating from the pulsar - these latter presumably aligned 
with its rotation axis and proper motion vector (Caraveo & Mignani 1999; Ng & Romani 
2007; Kaplan et al. 2008 and references therein). The "inner-ring", prominent in X-rays, 
is commonly accepted as being the termination shock produced by the relativistic wind of 
particles accelerated by the pulsar. Many of the optical and X-ray features brighten and 
fade and/or move over weeks or months (e.g., Hester et al. 1995; Hester et al. 2002). 

The quiescent or average spectral energy distribution (SED) of the Crab Nebula has a 
characteristic two- humped form (see, e.g., Figure El Atoyan & Aharonian 1996; Bucciantini, 
Arons, & Amato 2011 and references therein). The synchrotron spectrum extends from 
« 30 MHz to « 1.2 x 10 22 Hz (500 MeV). Most of the power is radiated by « TeV electrons 
in the near UV ~ 10 eV with an associated luminosity of ~ 1.3 x 10 38 ergs/s (Hester 2008). 
This roughly matches the loss of rotational energy by the pulsar, which releases its energy 
electromagnetically, generating a current ~ 200 TA and inducing an electro-motive force 
(EMF) fa 50 PV. However, the nebula is currently varying on a few-year timescale (Wilson- 
Hodge et al. 2011). At higher energies, Compton scattering has a luminosity of ~ 10 36 erg/s, 
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peaking around 60 GeV (Albert et al. 2008) and measured up to ~ 80 TeV (e.g., Aharonian 
et al. 2004; Abdo et al. 2010). 

Since 2007, the AGILE and Fermi satellites have detected several 7-ray flares from the 
Crab Nebula (Tavani et al. 2011; Abdo et al. 2011; Striani et al. 2011a; Buehler et al. 2012) 
in the 0.1 — 1 GeV range. The most dramatic flares exhibit variability on timescales as short 
as a few hours, although it is unclear whether they are distinct events or just the largest 
variations from a stationary power spectrum of fluctuations. Prior to the 2011-April event, 
the only Crab 7-ray flare covered by a multi-wavelength observing program was the 2010- 
September flare, which triggered observations in radio, optical (using both ground-based 
telescopes and HST), and X-ray bands. Despite the 7-ray brightness of the flares, there has 
been no evidence for correlated variations in radio (Lobanov, Horns, & Muxlow 2011; this 
paper), near infrared (Kanbach, et al. 2010, this paper), optical (Caraveo et al. 2010), or 
X-ray bands (Evangelista et al. 2010; Shaposhnikov et al. 2010; Tennant et al. 2010; Ferrigno 
et al. 2010; Horns et al. 2010; Cusumano et al. 2011; Tennant et al. 2011; Tavani et al. 2011; 
Striani et al. 2011b; this paper). 

Here we focus on the Fermi-LAT results for the 2011-April flare (Buehler et al. 2012), 
which allow us to assess the source behavior in detail. The source doubled its 7-ray flux 
within eight hours and reached a peak flux 30-times its average. The isotropic luminosity 
increased to w 2 x 10 37 erg/s in « 10 hr and the spectrum peaked at ~ 400 MeV. Tabled] 
gives the 7-ray powerlaw photon spectral index, the integrated photon flux above 100 MeV, 
and the photon spectral flux at 100 MeV, as measured during the 10-ks time intervals when 
X-ray data ($2$ were taken. 

Notification as to the level of flaring prompted us to trigger pre-approved Target of 
Opportunity observations with Chandra and with the NRACo Karl G. Jansky Very Large 
Array (VLA). We were also fortunate to obtain a Keck image in the near infrared, albeit 
not under ideal conditions. Figure [1] shows the Fermi-LAT 7-ray counting rate as a function 
of time and also indicates the times of the Chandra, Keck (§|3]) and some of the VLA (§4|) 
observations. 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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2. X-ray Observations and Data Analysis 

With the back-illuminated ACIS S3 CCD on the Chandra X-ray Observatory approx- 
imately centered on the Crab pulsar, we obtained five observations (Table [1]) during and 
somewhat after the 2011-April 7-ray flare. For these observations, the spacecraft dithered 
with an amplitude set to 1". Although standard processing typically produces an aspect 
solution better than 0.5", even this small uncertainty can introduce noticeable shifts when 
comparing different data sets. Thus, we re-registered images for our analysis using the 
read-out streak and the pulsar as guides. As each of these 5 images was placed at approx- 
imately the same CCD location, spatial non-uniformity in the ACIS response (e.g., due to 
contamination) does not introduce spurious temporal variability. 

Owing to the Crab's high flux, the ACIS observations employed a special mode with 

0. 2-s frame time, which limits the CCD read-out to a 300 x 300 ACIS-pixel (w 150" x 150") 
subarray. Although each observation lasted about 10 ks, telemetry saturation reduced the 
effective integration time to approximately 1200 s per observation. Despite the short frame 
time of the special ACIS mode, regions of high surface brightness suffer somewhat from 
pile-up effects. We consider only data in the range 0.5-8.0 keV because of severe interstellar 
absorption at low energies and declining flux at high energies. Using these data, we then 
search for X-ray variations approximately contemporaneous with the 2011-April 7-ray flare. 

2.1. X-ray Image Analysis 

Figure |2] shows an image of the number of counts per ACIS pixel, summed over the 
5 observations. For each observation, we re-binned a 120 x 120 ACIS-pixel image centered 
on the pulsar into a 60 x 60 array of 2 x 2 ACIS pixels. Each of these I = 3600 "analysis 
pixels" is sufficiently large (about 1 square arcsec) to enclose most of the Chandra point 
spread function anywhere in the field of view. Note that we also performed an analysis using 
a circular bin of radius 1 ACIS pixel on an oversampled grid of cadence 0.1 x ACIS pixel - 

1. e., a spherical top hat smoothing of the events. As each method gave similar results, we 
here report the results for the "analysis pixel" binning, for which each pixel is statistically 
independent. 

For each analysis pixel i, we calculate the mean count rate averaged over the J = 5 
observations, weighted^ by the respective (counting-rate) statistical error er^ for each analysis 
pixel and observation. For evaluating statistical significance of temporal variations over the 
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J = 5 observations, we computed Xi along with a derived significance measure Si. 

For purposes of discussion, we also compute the appropriately weig htecfl statistical error 
(Tj and sample standard deviation Sj for each pixel i. As properly weighted, Xi — (J~ 1) s ?/°f • 
While we have rigorously calculated r», xf, CT i> anc l s % f° r eacn pixel z using appropriate 
weightings, we note that the weightings are nearly uniform as the effective duration of the 
each of the J = 5 observations was nearly the same — about 1200 s. 



2.2. Variability of the X-ray images 

The (counting) statistical error <j{ is the primary noise term and thus governs the sensi- 
tivity for detecting temporal variations at the analysis-pixel ( square- arcsec) scale. Figure [3] 
shows the image and the corresponding histogram of the distribution of Oi, which ranges 
from 0.0025 to 0.024 ct/s per analysis pixel. Based upon the x 2 probability distribution and 
the number of "tries" (J = 3600 independent analysis pixels), a 99%-confidence detection 
would require a X?99% > 31.2 on (J — 1) = 4 degrees of freedom. This corresponds to a 
sample standard deviation > 2.80 o"j, which ranges from 0.0071 to 0.068 ct/s over the 

field. 

We do not here display the analogous image and histogram for the sample standard 
deviation Si, which ranges from 0.0014 to 0.048 ct/s per analysis pixel. Instead, Figure H] 
shows the image and histogram of the distribution of a calculated ( §2.ip significance measure 
Si, related to xl — 1) s i/°f- The statistically most significant variation has xt — 23.5 on 
v = (J— I) = 4 degrees of freedom giving Si = 6.9. Such a fluctuation is expected statistically 
in at least 1 of 3600 pixels in 31% of realizations. Table[2]gives the sample standard deviation 
and 99%-confidence upper limit to the count-rate variation, for the analysis pixel with the 
statistically most significant X-ray variation. While we detect no variations statistically 
significant at 99%confidence, it is curious that the 3 most significant variations occur at 
locations on the inner ring. 

Note that if a feature, such as one of the knots, doesn't change in intensity but moves 
from one analysis pixel to another, then our x\ test would detect this as a variation. We know 
that features in the inner ring of the Nebula do move and expect to detect some variability 
due to this motion. However, as our 5 observations span only 14 days and 1" corresponds to 
11.5 light days, only relativistic motion would be detectable. Other effects, such as changes 



3 X? = E/=i{(r« " n)74> and S * = bd - where n = ( J - 1). 

4 of = J/E/=i{i/4> and s " = iAEn{(ni -^)74}/£^i{i/4} = A a * 
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in the roll angle of the read-out streak, can also lead to spurious variability. Indeed, this 
may play a role for the analysis pixel with the most significant variation, which lies adjacent 
to the average read-out streak (Figure [2]). 

2.3. Limits to the X-ray flux 

Thus far, we have described the X-ray data for each analysis pixel in units of ACIS 
count rate. Neglecting for the moment pile-up effects, the photon spectral flux (or other 
related radiation quantity) is proportional to the count rate for an assumed spectral shape. 
Consequently, any change in count rate corresponds to a proportionate change in photon 
spectral flux (for an assumed spectral shape). Using the Chandra PIMM £@ for the ACIS-S 
detector and an absorption column A% = 3.1 x 10 21 cm -2 , we determine (ignoring pile-up) 
this constant of proportionality for an X-ray power-law photon index T x = |, 1, and 2: At 
E x = l keV, N E (E x )/r = 0.99, 1.26, and 2.46 xl(T 3 ph/(cm 2 s keV) per ct/s, respectively. 
Correcting for pile-up has little effect in low-count-rate regions, but would raise these flux 
upper limits by ~ 10% or so for high-count-rate regions. 

Table [2] calculates the photon spectral flux N E (E X ), the energy spectral flux F E (E X ), and 
the indicative (isotropic) luminosity EL E (E X ) = 4ttD 2 EF e (E x ) at D = 2 kpc, corresponding 
to the sample standard deviation and 99 %-confidence upper limit for the count-rate variation 
in the analysis pixel with the most significant X-ray variation. Figure |5] displays an image of 
the energy spectral flux F E (E X ) at E x = 1 keV for = 1, based upon the sample standard 
deviation Si of the count rate in each analysis pixel. 

2.4. Constraints on the X-ray to 7-ray Spectral Index 

We now compare the X-ray data with the 7-ray data to quantify the implications of our 
lack of detection of time variations in the X-ray data. Our approach compares a variability 
measure for the X-ray (1-keV) photon spectral flux AN E (E X ) in each analysis pixel with the 
analogous variability measure for the 7-ray (100-MeV) photon spectral flux AN E (E- r ). In 
particular, we calculate the sample standard deviation of the 7-ray spectral flux at 100 MeV, 
using power-law fits to the 5 Fermi-LAT measurements that were simultaneous with the 
5 Chandra observations. Table [1] lists the 5 Chandra ObsIDs, their dates, along with 
the 7-ray photon index T 7 , integrated photon flux iV(>100 MeV), and photon spectral flux 



5 http://asc. harvard.edu/toolkit/pimms.jsp 
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Ne (100 MeV). For the 5 Fermi-L AT observations, the mean and sample standard deviation 
of the photon spectral flux at 100 MeV are 1.21 x 10~ 10 and 5.77 x 10" 11 ph/(cm 2 s keV), 
respectively. 

Based upon the sample standard deviation (sj) of photon spectral flux at E x = 1 keV 
for each X-ray analysis pixel and the measured standard deviation (5.77 x 10 -11 ph/(cm 2 s 
keV)) at Ery = 100 MeV, we constrain the effective X-ray to 7-ray photon index of the flaring 
component: Y X1 = — \og[AN E (E^) / AN E (E X )}/ \og[E 7 / E x }. Figure |6] shows the image and 
corresponding histogram of the distribution of upper limits to r x7 based upon the sample 
standard deviation of the X-ray measurements and assuming T x — 1. 

In that the 7-ray variations are statistically significant and the X-ray variations are 
not, we compute 99%-confidence upper limits to T X1 (Table [2] last row). Note that the 
upper limits to T X1 are marginally consistent with the low-energy extrapolation of the 7-ray 
spectrum (T 7 = 1.27 ± 0.12) of the flaring component (Buehler et al. 2012). 

2.5. Variability within an X-ray Image 

In sections I2.1H2.2} the search for variability focused on sensitivity to flux changes 
amongst the five pointings with a minimum cadence of 0.6 days. Here, we search for vari- 
ability on shorter time-scales — namely within each pointing. As described in Section [2j 
the ACIS S3 CCD was read at most roughly 6000 times in each pointing due to telemetry 
saturation (deadtime). In this study, rather than using the 2x2 ACIS analysis pixel, we 
employed a circular search bin with a radius of 1 ACIS-pixel (0.49") on a grid with 0.1 ACIS 
pixel spacing. Note that this oversampling implies that the results of the test in adjacent 
pixels are not statistically independent. (We also analyzed these data using the statistically 
independent analysis pixels of §2.11 with similar results as below.) 

Using the frame number of each detected photon, we derive the empirical cumulative 
distribution function (ECDF) of the frames with a photon arriving in the analysis pixels 
of the CCD. This ECDF is then compared with the corresponding ECDF of the exposure 
given by the sequence of frames actually read. Finally, we compare the two ECDFs using 
a Kolmogorov-Smirnov test, resulting in a probability estimate Q that the two ECDFs 
are derived from the same parent distribution. A low value of Q would indicate possible 
variability. The results of the test for the five pointings were very similar. The smallest 
value, Q m i n = 6.7 x 10~ 7 , was obtained in observation 13151. Note that selecting this point 
represents a tuning bias, as the noise in neighboring points is highly correlated due to the 
oversampling. The probability of finding at least one pixel with Q m i n considering that there 
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are 120 x 120/(7r x l 2 ) x 5 statistically independent trials is 0.015, which we regard as a 
lower limit due to the tuning bias. A 0.015 probability is tantalizing but not compellingly 
significant: Hence, we do not claim detection of short-time-scale variability. The fact that 
the location of the point with minimum Q is very close to the pulsar, a region in which pileup 
plays a strong role in blotting out the image, bolsters our somewhat conservative conclusion. 

3. Near-Infrared Image of the Inner Knot 

The extreme saturation of the pulsar in the X-ray images means that we cannot easily 
study the central 2" in X-rays. However, this region does contain a nebular structure of 
particular interest: the "inner knot" whose peak is 0.65" southeast of the pulsar at position 
angle 118° East from North (Hester 2008). This structure, an oval shape extending 0.75", 
is well measured in HST and ground-based near-IR images. Given its relatively red spectrum 
(energy spectral index a v = —1.3 ±0.1 versus a v = 0.27 ±0.03 for the pulsar; Sandberg and 
Sollerman 2009), it is one of the near-IR brightest structures in the Nebula. Sandberg and 
Sollerman (2009) note that the knot varies by a factor of 2; we confirm typical variability of 
20 — 30% in archival HST images. Komissarov & Lyutikov (2011) have proposed that this 
structure represents radiation from an oblique termination shock in the pulsar wind nebula. 
In this picture, the Earth line-of-sight is tangent to the flow at the inner knot position, 
and thus the intensity experiences substantial Doppler boosting for synchrotron emission in 
the mildly relativistic post-shock flow. Indeed, in relativistic MHD simulations they find 
that this bright spot is highly variable and can dominate the 7-ray synchrotron emission. 
Alternatively, the knot could be a time varying standing shock in the polar jet flow itself, a 
flow known to be highly variable from HST imaging (Hester 1995, 2002, 2008). 

It is thus of interest to check the status of the knot during the 2011- April 7-ray flare. 
Unlike the sequence of multiwavelength observations performed after the 2010 September 
flare, it was impossible to trigger an allocated HST Target of Opportunity observation owing 
to solar constraints in April. Happily we were able to obtain a Keck Near Infrared Camera 
(NIRC2) K' exposure (Figure [7] left image) on MJD 55667.250, almost precisely at the peak 
of the 7-ray flux and 2.5 h before the ACIS image ObsID 13152 (Figure [T]). Unhappily, the 
observations occurred during twilight and only one 20 x 4-s integration without dithering 
was obtained. Under these conditions the adaptic-optics (AO) loop did not close, leaving an 
undithered image with native 0.46" full width at half maximum (FWHM) seeing. This frame 
was dark subtracted and an approximate background was removed using an immediately 
subsequent image. Despite the modest image quality, the inner knot was well detected. 
After subtracting the pulsar with a scaled image of the comparably bright companion star 
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4" northeast, we measured the knot flux and position. We find a magnitude K' = 15.60±0.03 
and an offset 0.64±0.04" from the pulsar. For comparison we measured a high-quality NIRC2 
K' image (Figure [7] right image) obtained 2005 Nov 10. Here the knot is K' = 15.94 ± 0.02 
at offset 0.58 ± 0.02". We also note that Sandberg and Sollerman (2009) measured K s = 
15.80 ± 0.03 on 2003 Oct 18. We conclude that the knot was in a relatively bright state 
during the flare (~ 35% brighter than in 2005), but well within the normal range of flux 
(and position) variation. Thus, there is no dramatic change in the inner knot in the near-IR 
band. We use the amplitude of the measured variation as an upper limit to any variation in 
the inner knot associated with the 7-ray flare (Figure [8]). 



4. Radio Observations 

On 2011 April 14, we triggered a prompt radio follow-up program with the VLA. The 
VLA observations occurred in 8 epochs starting April 15 and ending July 10. These ob- 
servations detected the pulsar at 2 epochs, but found no other point sources in the field. 
Observations were predominantly in the range 4 — 8 GHz, with additional observations at 
1.4 GHz (not reported; see below) and at 22 GHz (for some later epochs). Unless otherwise 
noted, all observations used two sub-bands, each with a 128-MHz bandwidth. In each run, 
observations of the target were bracketed with scans of a phase calibrator (J0559+2353, ex- 
cept where noted) and a flux calibrator (3C 147). The fields of view are limited to the primary 
beam response of the antennas, with full width at half power of 9'/^ with 1/5 = v/(5 GHz). 
Table [3] provides a summary of the observational parameters and results. 

Our initial observations were obtained through a Fermi guest-investigator cycle-3 pro- 
gram (S3184) approved for four 1-hour runs in the L-band and C-band (~ 1.4 and 5 GHz, 
respectively). The VLA was in its B-array configuration during these observations, resulting 
in images with angular resolution « 1"/ z/5. We found that the L-band data for the Crab were 
highly confused due to the brightness and complexity of the steep-spectrum nebular emission 
in the first (April 15) and second (April 19) epochs. Consequently, we modified our strategy 
for subsequent observations. After the first epoch, we split the C-band observations into 
two widely spaced side-bands centered at 4.2 and 7.8 GHz, aiming better to constrain the 
spectrum of any detected source. We also began scheduling observations only at frequencies 
greater than 4 GHz after the second epoch. In these B-array data, our point-source limits 
at the lower frequency are ~ 3 — 4x larger than at the higher frequency, again due to the 
Crab Nebula's steep-spectrum radio emission. 

After non-detection of any significant radio point-source emission down to w 1-7 mJy 
(3-sigma) sensitivities in the initial three VLA observations 1-7 days after the 7-ray peak 
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(Hays et al. 2011), we purposely delayed the fourth observation until 9 days after the previous 
observation to probe longer time scales. Only in this last observation (April 30) did we obtain 
a significant point-source detection, which was coincident with the Crab pulsar position, but 
only in the upper side-band centered at 7.8 GHz. The detection is a factor of 5 greater than 
the 3-sigma limit from 9 days prior. The source was not detected at the lower frequency 
side-band (4.2 GHz) with a limit indicating a source with a flat radio spectrum. 

Following the point-source detection on April 30, we became aware of VLA TEST 
observations of the Crab obtained on April 22 (program TDEM0007, PI: D. Frail). These 
data were obtained with wide bandwidth (16 x 128 MHz wide sidebands), so were more 
sensitive than those from our observing runs. The flux densities were scaled to 3C 147; 
J0534+1927 was utilized for phase calibration. The flat-spectrum radio source coincident 
with the Crab pulsar detected in our Apr-30 observation was confirmed in the Apr-22 data 
in two bands, but with a much lower (10 x) flux. Also, the source spectrum was rather steep, 
with an energy spectral index a = 2.21 ± 0.34 (S u oc v~ a ) between 5 GHz and 8.6 GHz. 

Following the radio detections of the pulsar, we requested further VLA monitoring of 
the Crab through Director's Discretionary Time (program 11A-268 = AC1052). In addition 
to the C-band observations, we obtained exposures in the K band (centered at 22.396 and 
22.254 GHz) aiming to constrain further the spectrum of any detected radio source. Through 
this program, we obtained 2-hour runs on May 12/13 (while the VLA was in its hybrid BnA 
array) and on July 10/11 (in A array), and an additional 1-hour run on May 28, using one 
of the early (April 19, from program S3 184) frequency setups. An angular resolution of 
~(0.3"/z/5) is typically achieved in A-array VLA observations. With the higher resolution, 
we obtained systematically 4x lower flux limits than in the lower resolution B- array data, 
presumably due to lesser contribution from the extended nebular emission. In none of these 
later epoch follow-up observations, did we detect a point source, to typical limits of 1 — 2 
mJy at each of the three frequencies (see Table 



4.1. Discussion of the Radio Data 

Previously, it was argued that the 7-ray flaring possibly originates in a knot 5.7" east 
of the pulsar (Tavani et al. 2011). Indeed, this knot is the site of the most significant X-ray 
variability we observe (§ 12. 2p during the 2011-April flaring episode. Variable radio emission 
was detected around the time of the previous Crab 7-ray flaring episode (Lobanov, Horns, & 
Muxlow 2011) with fainter flux densities than achieved in our VLA observations. However, 
we found no significant radio point-source counterpart to this knot in any of our 8 epoch VLA 
observations following the 2011-April 7-ray flare. Rather, we detected a variable continuum 
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radio source with the VLA, coincident within 0.2" of the Crab pulsar position in 2 of 8 
epochs. However, the flux level and cadence of the radio detections is consistent with previous 
observations by Moffett & Hankins (1996), who detected the pulsar 20 — 40% of the times 
they observed. Moreover, the dates of the radio point-source detections coincident with the 
pulsar do not coincide with any feature in the 7-ray lightcurve (Figured]), having occurred 
8-16 days after the brightest 7-ray peak. Consequently, our VLA follow-up observations 
provide no conclusive evidence for the site of the 7-ray flares. 

5. Discussion 

Here we discuss possible explanations for the absence in non-7-ray bands of variability 
that is obviously correlated with the 7-ray flare. In addition we present a conceptual model 
for the production of the 7-ray flares. 

5.1. 7-ray Emission 

As was recognized immediately, the SED of 7-ray flares peak near a characteristic energy 
about 5 times the energy a _1 m e c 2 ~ 70 MeV, which is identified with radiation-reaction- 
limited synchrotron (magneto-bremsstrahlung) emission (e.g., Landau & Lifshitz 1959). Sub- 
jected to comparable parallel and perpendicular electromagnetic acceleration, an electron 
radiates at this energy, independent of the strength of the acceleration. An electron emitting 
synchrotron radiation in a magnetostatic field with peak emission at several 100 MeV would 
cool in turning through « 0.2 radian and would thus require a parallel electric field E w 5cB 
to compensate the radiative loss. 

The Crab pulsar releases energy in an essentially electromagnetic form. Poynting flux 
flows radially outward from the the pulsar through the light cylinder at cPj (2ir) w 1500 km 
and into an outflowing wind, where at least some of the electromagnetic-energy flux may 
transform into a plasma-energy flux. How, where, and to what extent this happens has long 
been a matter of debate (e.g., Arons 2010; Kirk et al. 2009). Furthermore, the electromag- 
netic component has a DC toroidal part with an associated quadrupolar current distribution, 
and an AC, "striped" part containing current sheets separated by |cP ~ 5000 km. The trans- 
formation from electromagnetic to plasma energy might be non-dissipative — through the 
action of a Lorentz force (e.g., Bogovalov 1997; Bogovalov 2001) — or dissipative — through 
particle heating and acceleration (e.g., Coroniti 1990; Lyubarsky & Kirk 2001; Sironi & 
Spitkovsky 2011). However, it must occur somewhere as magnetic flux would otherwise 
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accumulate in the nebula, ultimately reacting back on the pulsar. Some of this transfor- 
mation from electromagnetic to plasma energy may occur at a shock (Petri & Lyubarsky 
2007; Sironi & Spitkovsky 2011) with radius ~ 10 17 cm, where the wind momentum flux 
balances the ambient nebular pressure (Rees & Gunn 1974; Kennel & Coroniti 1984). It has 
also been proposed that the toroidal field loops contract to form an axial pinch (identified 
with the X-ray jet) and reconnect at an equatorial current sheet (the torus) (Komissarov & 
Lyubarsky 2003; Del Zanna, Amato, & Bucciantini 2004; Camus et al. 2009). 

In many respects the pulsar is a current generator. The supersonic wind contains out- 
flowing fluxes of electrons and positrons. (Any ions that are present behave like positrons 
of similar rigidity but do not radiate.) The difference in their fluxes determine the current 
density. This current may concentrate into sheets and filaments, where strong dissipation 
can occur — as happens in heliospheric and laboratory plasmas (Gosling et al. 2005; Sui & 
Holman 2003; Sergeev et al. 1993). In particular, the inner wind, the shock, the jet, and 
the torus are all natural sites of rapid dissipation and 7-ray emission. If we consider this 
dissipation more generally under electromagnetic conditions, a current / may be associated 
with a potential difference V ~ IZ where Z = fi c = 377 Q is the impedance of free space 
and we drop model-dependent constants of order unity. (This result can be anticipated on 
the basis of dimensional analysis or exhibited in particular simple cases.) The maximum 
energy to which an electron or positron can be accelerated is 7max^eC 2 ~ eV ~ eIZ and the 
expected power is then L pa IV ~ I 2 Z . On this basis, a spectrum of currents extending up 
to pa 30 TA should suffice to account for the 7-ray variations. 

However, currents do not automatically dissipate as just described. Large electric fields 
are normally discharged in a few plasma periods. The best way to create them here is tran- 
siently over a few Larmor periods and radii. This, in turn, requires local charge separation 
of the plasma in the emission site. To be more precise, the density n = n_ + n + of electrons 
plus positrons/ions will combine to create a local current density j ~ nec. For example, in 
the case of a pinch, the gradient, polarization and curvature drifts automatically produce 
the axial current. However, a supporting local electric field of strength E pa cB also requires 
that |n_ — ri+l pa n. 

Put another way, the charge and current are mostly in an emission site that is a few 
Larmor radii in size and survives for a few Larmor periods of the 7-ray emitting particles. 
When this happens, the "Ohmic" dissipation is radiative, not collisional as is normally the 
case. For this to occur, the particles must be sufficiently energetic to radiate efficiently. This 
requires that most of the particles are concentrated in an emission site that is as small as 
pa 7 r 3 nax r e pa (eIZ /m e c 2 ) 3 r e pa 10 16 cm. The key point is that there should be extensive and 
sustained radiation-reaction-limited emission at the peak 7-ray-flare energy, even though 
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the particle energy and magnetic field might be changing. Note that when this condition 
is unsatisfied, efficient particle acceleration to lower energy should still result: Most of the 
magnetic dissipation and particle acceleration in the nebula might occur in this fashion. 
Detailed modeling is necessary to determine whether or not such a scheme can reproduce 
the powerful, narrow-band 7-ray variation that is observed (Uzdensky et al. 2011; Cerutti et 
al. 2012; Cerutti, Uzdensky, & Begelman 2012; Lyutikov, Balsara, & Matthews 2012; Bykov 
et al. 2012; Sturrock & Aschwanden 2012; Blandford & Yuan 2012) and to see if unstable 
magnetized plasmas, carrying large currents, evolve to satisfy these conditions. 

5.2. Associated Emission 

Whether we interpret the 7-rays as coming from radiation-reaction-limited synchrotron 
emission or simply extrapolate the observed 7-ray spectra to lower energy, it should not be 
surprising that direct, associated emission has not yet been observed in the X-ray, optical, or 
radio bands: In these bands the contrast with the steady emission is too small to be easily 
noticed. However, the indirect effects could be larger and detectable. For example, the large 
2011- April flare produced a radiant energy of 6 x 10 40 ergs if isotropic, equivalent to the 
energy contained within a region of size ^2x 10 16 cm subtending an angle m 0.3 arcsec. It 
seems unlikely that the dynamical aftermath of a major flare would not alter the ambient 
emission — either through compression or rarefaction that would cause the magnetic field 
strength and the electron distribution function to change significantly. The associated surface 
brightness change should be several percent, assuming a total emission region of size ~ 
0.3 arcsec, consistent with our upper limits. Even if future observations fail to exhibit 
associated emission, they may still rule out specific detailed mechanisms in local sites. 

Understanding the emission mechanism could have a significance beyond pulsar wind 
nebulae. In particular, it could provide a clue to the surprisingly rapidly variable emission 
seen in relativistic jets in the radio, optical, X-ray, and TeV bands. If so, the Crab Nebula 
would once again be the source of fresh and important astrophysical insight. 

6. Summary 

Using the Chandra, Keck, and VLA Observatories, we acquired X-ray, near-IR and radio 
images of the Crab Nebula, contemporaneous with the 2011- April 7-ray flare. We searched 
for variability in the X-ray data over two time-scale ranges: First we tested for pointing-to- 
pointing variations amongst the 5 pointings, each with an effective exposure time ~ 1200 s 
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and a minimum separation of 0.6 days. Second we tested for variations within each of the 
5 observations. In neither case did we detect statistically significant X-ray variations; thus 
we can set only upper limits to any X-ray variations associated with the 7-ray flare. As the 
Chandra ACIS images suffer severe pile-up near the Crab pulsar, our search for variability 
in the X-ray images was not sensitive to variations within the central ~ 1.5" or so. 

Comparing the upper limits to X-ray variations with the Fermi-LAT-measured 7-ray 
variations, we set upper limits at 99%-confidence to the effective X-ray-7-ray photon power- 
law index T XJ < 1.20 to < 1.27, dependent upon assumptions about the X-ray index T x . 
As Fermi-LAT measures a 7-ray index T 7 = 1.27 ± 0.12 for the flaring component, it is 
statistically possible that the flaring component's spectrum extends as a simple power-law 
from 7-rays to X-rays. Further, we note that our upper limit to is consistent with 
transparent synchrotron emission, whose photon index must be > |. 

Comparison of two Keck near-IR observations found that the inner knot (~ 0.65" from 
the pulsar) was somewhat brighter than average during the 7-ray flare, but well within the 
normal range of brightness fluctuations typically observed. We used the measured (~ 35%) 
change in the near-IR flux from this knot as an upper limit to near-IR variations associated 
with the 7-ray flare. We also performed a number of VLA observations searching for a point 
source appearing either at an unusual location and/or contemporaneous with the 7-ray flare. 
Other than the pulsar itself, no such source was detected. 

Figure [H] shows the spectral energy distribution (SED) of the Crab Nebula over the 
observed electromagnetic spectrum. The plot also shows the SED of the 2011-April 7-ray 
flare and the various limits determined here on variable radio, near-infrared, and X-ray 
emission possibly associated with the 7-ray flare. 

Finally, we reviewed and discussed potential implications of 7-ray flares and theoretical 
issues to be addressed. We concluded that, apart from lower-energy emission directly asso- 
ciated with the 7-ray flare itself, the dynamical aftermath of a major flare could alter the 
ambient emission — e.g., through compression of the magnetic field. The associated surface 
brightness change would likely be only several percent, assuming a total emission region of 
size ~ 0.3", consistent with our upper limits. 

Although no "smoking gun" has been identified, one should be encouraged that we have 
identified a number of regions in the X-ray images that are possible candidates. We have 
also established further Target of Opportunity observations with Chandra and HST that 
will be triggered at the onset of the next 7-ray flare. The X-ray observations will also probe 
the region very close to the pulsar using the Chandra High- Resolution Camera (HRC). 
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Table 1 : Time-ordered list of Chandra observations and Fermi-L AT 7-ray measurements' 1 . 



ObsID 


Date b 


r c 


N(> 100 MeV) d 


N E (\m MeV) 


e 


13150 


4851.039 


2.42 ± 0.08 


(1.27 ±0.08) x 10" 5 


(1.80 ±0.15) x 10" 


10 


13151 


4851.667 


2.25 ±0.15 


(8.05 ±0.15) x 10~ 6 


(1.01 ±0.12) x 10- 


10 


13152 


4853.423 


2.20 ±0.06 


(1.54 ±0.06) x 10~ 5 


(1.84 ±0.12) x 10- 


10 


13153 


4859.032 


2.27 ±0.21 


(4.68 ±0.21) x 10~ 6 


(6.0 ± 1.0) x 10" 


11 


13154 


4865.335 


2.76 ± 0.40 


(4.50 ± 0.40) x 10" 6 


(7.9 ± 1.9) x 10- 


11 



a Analyzed following all the procedures in Buehler et al. (2012) 

b Days after MJD 50814 to the middle of the observation, which is 10-ks long. 

c 7-ray powerlaw number index. The error is the maximum of the two-sided uncertainty 

d Photon integrated flux [ph/(cm 2 s)] above 100 MeV. The error is the maximum of the 

two-sided uncertainty 

e Photon spectral flux [ph/(cm 2 s keV)] at 100 MeV 



Table 2: X-ray results at 1 keV for the analysis pixel with the most significant variation. 



Quantity 


Unit 


Sample stdev (s) 


99%-upper 


limit 


Rate 


ct/s 


0.0480 


0.0554 




r x 




1 1 2 


2 1 

3 


2 


N E 


10- 4 ph/(cm 2 s keV) 


0.48 0.61 1.18 


0.55 0.70 


1.36 


F E 


1Q-13 er g/( cm 2 s keV) 


0.76 0.97 1.89 


0.88 1.12 


2.18 


ELe 


10 32 erg/s 


0.37 0.47 0.91 


0.42 0.54 


1.05 


r 

j- £7 




1.18 1.20 1.26 


1.20 1.22 


1.27 
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Table 3: VLA point-source detections and limits in the Crab field 



Epoch 


Date a 


Array 


Frequency 


Detection or Limit c 


(in 2011) 






(GHz) 


(mJy) 


Apr 15 


4852.980 


B 


4.959 


< 7.1 


Apr 19 


4856.006 


B 


4.195 


< 5.7 


Apr 19 


4856.006 


B 


7.795 


< 1.7 


Apr 21 


4858.011 


B 


4.195 


< 4.3 


Apr 21 


4858.011 


B 


7.795 


< 1.0 


Apr 22 


4859.919 


B 


4.910 


1.78 ± 0.086 


Apr 22 


4859.919 


B 


8.566 


0.52 ± 0.094 


Apr 30 


4867.901 


B 


4.195 


< 7.9 


Apr 30 


4867.901 


B 


7.762 b 


5.55 ± 0.86 


Apr 30 


4867.901 


B 


7.827 6 


4.59 ± 0.96 


May 12 


4879.985 


BnA 


4.195 


<2.1 


May 12 


4879.985 


BnA 


7.795 


<0.8 


May 12 


4879.985 


BnA 


22.46 


<0.4 


May 28 


4895.925 


BnA 


4.195 


<2.4 


May 28 


4895.925 


BnA 


7.795 


<1.2 


Jul 10 


4938.791 


A 


4.195 


<0.7 


Jul 10 


4938.791 


A 


7.795 


<1.1 


Jul 10 


4938.791 


A 


22.46 


<1.7 



Notes - 

a Days after MJD 50814 to the middle of the on-target exposures. 

b For the 7.8 GHz detection on Apr 30, the data were further split into two sub-bands. 

c Limits list 3-sigma uncertainties. Detections list 1-sigma uncertainties; corresponding 

point-source limits during these 2 epochs should take 3 times the 1-sigma values. 
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Fig. 1— Fermi-LAT photon flux (10~ 7 ph/(cm 2 s)) above 100 MeV during the 2011-April 
flare as a function of time. Displayed data extend beyond the time span shown in Buehler 
et al. (2012) but follow the same data processing as described there: Data are adaptively 
binned with a 20-minute average bin duration. The full-range vertical lines denote times of 
the 5 (~ 1200-s) Chandra observations (black), of the Keck observation (blue, fourth from 
left), and of the first 5 VLA observations (green). 
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Fig. 2. — Summed image for the 5 Chandra ACIS observations occurring near the 2011-April 
7-ray flare, at the native CCD resolution. The color bar gives summed counts per ACIS pixel 
over a total effective exposure of about 6 ks. North is up and the pulsar is at (0,0) in the 
displayed ACIS-pixel coordinates. The nearly horizontal read-out streak through the pulsar's 
location is the trailed (out-of-time) image, resulting from exposure of each CCD pixel as the 
image is read out at 40 /is per row. As the 5 observations occurred at slightly different roll 
angles, the read-out streak is slightly blurred azimuthally. The X symbols mark locations of 
the 3 statistically most significant variations (Si > 6, §2.21) . the most significant lying to the 
east of the pulsar. 
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Fig. 3. — Left: Image of the statistical error <jj (ct/s) in the counting rate per analysis pixel 
(about 1 square arcsec) for the 5 Chandra ACIS observations. North is up and the pulsar 
is at (30,30) in the displayed analysis-pixel coordinates. Right: Histogram of number of 
occurrences of each value in the image to the left. 
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Fig. 4. — Left: Image of the significance measure Si = (xf — v^)l^f2vi of sample counting- 
rate variations amongst the 5 observations. North is up and the pulsar is at (30,30) in the 
displayed analysis-pixel coordinates. Right: Histogram of number of occurrences of each 
value in the image to the left. 
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Fig. 5. — Image of the energy spectral flux Fe{E x ), in erg/ (cm 2 s keV) at E x = 1 keV, based 
upon the sample standard deviation (sj) of the counting rate and assuming T x = 1. Note 
that the indicative energy flux F(E) = EFe{E) in erg/(cm 2 s) happens to have the same 
numerical value as Fe(E) at E = 1 keV. North is up and the pulsar is at (30,30) in the 
displayed analysis-pixel coordinates. 
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Fig. 6. — Left: Image of upper limits to the effective photon index T X1 between 1 keV and 
100 MeV, based upon the sample standard deviation (s») of the counting rate and assuming 
T x = 1. North is up and the pulsar is at (30,30) in the displayed analysis-pixel coordinates. 
Right: Histogram of number of occurrences of each value in the image to the left. 




Fig. 7. — Keck K' observations of the Crab, after subtraction of a field-star image from the 
pulsar position (marked with x). Left: (MJD-50814) = 4853.250 without adaptive optics 
(AO). Right: (MJD-50814) = 2870, with laser guide star AO. Residuals from imperfect 
subtraction of the point spread function are visible at the pulsar position; the 'inner knot' 
is the extended structure to the southeast. 
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Fig. 8. — Spectral energy distribution (SED) of archival data (purple, orange, cyan, light 
green, turquoise, light blue) compiled by Meyer, Horns, & Zechlin (2010). Power is scaled 
from flux assuming isotropic emission at 2 kpc. The Fermi-LAT data for the 2011-April flare 
component appear in dark blue (Buehler et al. 2012). The solid black and red curves are 
fits to the flare spectrum with a power-law extrapolation to lower energies of photon index 
T = 1.27 ± 0.12 (spectrum 7 in Buehler et al. 2012). The three downward blue arrows at 
log 10 v = 17.4 mark 99%-confidence upper limits to a variable X-ray component, in increasing 
ELe = vL u for r^, = |, 1, and 2 respectively. (NB: Values for ELe from Table [2] are 
multiplied by 2-\/2 to scale from a standard deviation to a peak-to- valley, for comparison with 
the plotted SED of the 7-ray flare.) The red downward arrow at log 10 v = 14.1 indicates an 
upper limit to infrared variability of the inner knot, determined from the difference between 
the two Keck images. Finally, the black downward arrow at log 10 v = 9.7 gives an upper 
limit to 5-GHz radio variability, based upon the April- 15 VLA measurement. 
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